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Gamma Ray Bursts in the HAWC Era 


Peter Meszaros 1,3 , Katsuaki Asano 2 , Kohta Murase 1 , Derek Fox 1 , He Gao 1 , Nicholas Senno 1 


Abstract Gamma-Ray Bursts are the most energetic explosions in the Universe, and are among the most promis¬ 
ing for detecting multiple non-electromagnetic signals, including cosmic rays, high energy neutrinos and grav¬ 
itational waves. The multi-GeV to TeV gamma-ray range of GRB could have significant contributions from 
hadronic interactions, mixed with more conventional leptonic contributions. This energy range is important for 
probing the source physics, including overall energetics, the shock parameters and the Lorentz factor. We discuss 
some of the latest observational and theoretical developments in the field. 


1 Introduction 

Gamma ray bursts have continued to astonish and puzzle us since the early 1970’s when they were first publicly 
discussed. After it became clear that they were extragalactic, and hence had to be the most powerful and con¬ 
centrated explosive events in the Universe, which furthermore concentrate most of their output at MeV or higher 
energies, they have become a favored laboratory for exploring the most extreme environments and the highest 
energy processes in Nature. As such, they are the targets for increasingly sophisticated observing instruments, 
the latest of which is HAWC. 


2 The basic GRB paradigm 


In the course of the past 20-25 years of investigations, it has emerged that gamma-ray bursts (GRBs) can be 
classified in two main groups, “long” GRBs (LGRBs) and “short” GRBs (SGRBs) (e.g. Vedrenne & Atteia 


p009| )). The former have spectral peaks on average below 0.5 MeV and y-ray light curves whose duration 
nominally exceeds 2 s, while the latter have spectral peaks on average above 0.5 MeV and durations which 
mostly are less than 2 s. The stellar progenitors of LGRBs, which are generally located in star-forming galaxies, 
appear to be massive stars (M ^ 28M 0 ), whose core collapses after they have gone through their nuclear fuel 
burning cycle, while the progenitors of SGRBs are most likely compact binaries (double neutrons stars or neutron 
star-black hole binaries) whose merger is driven by gravitational wave emission (e.g.|Woosley & Bloom](|2006|); 
Gehrels et al.| ( [2Q09l )). See Fig.[I| 


The mechanism by which these stellar cataclysmic events liberate such huge energies 10 54 erg) in such 
short times (down to milliseconds) is ultimately linked to the liberation of the gravitational energy involved in 
the compact stellar core or in the merging of two degenerate stars (in both cases in a region of order a few 
Schwarzschild radii corresponding to a few solar masses withing a dynamic time of order milliseconds), and its 
conversion into thermal energy plus a final compact object, which ends up as a black hole. This thermal energy, 
produced explosively, manifests itself initially as a hot fireball (virial temperatures in the MeV range) consisting 
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Fig. 1 Progenitors of long and short GRBs, after M. 


Ruffert and H.-Th. Jahnka (1996), duration distribution from 


Kouveliotou et al. 


of photons, electrons, positrons, nuclei and magnetic fields, as well as neutrinos and gravitational waves. The 
prompt neutrinos produced are predominantly thermal, in the 10-20 MeV range, similar to those expected in core 
collapse supernova, and amount to a fewx 10 53 erg of energy, while the gravitational wave emission is mainly in 
the 10 2 — 10 3 Hz range, amounting again to a fewx 10 53 erg for SGRBs, down to possibly much less for LGRBs, 
depending on the variability of the total gravitational quadrupole moment during the merger or collapse. Both 
these neutrinos and the gravitational waves escape in a matter of seconds or less - the progenitor star is essentially 
transparent to them. 

The very high optical depth, however, traps the fireball photons, e ± , nucleons and magnetic fields, and this 
trapped thermal energy leads to an accelerating expansion of the fireball, converting thermal into bulk kinetic 
energy. For an initial ratio of energy to rest mass 17 = Eq/Mqc 2 1 the ultimately resulting bulk Lorentz factor 
eventually can become very large, F ^ rj . The observation of GeV photons would imply thermal spectra below 


the pair threshold ~ 0.5 MeV, unless the expansion Lorentz factor is very large F ~ 10 2 — 10 3 (e.g. Piran (2004); 


Meszaros (2006); Gehrels et al. (2009)), which in turn requires very low baryon loads in the fireball. The ob¬ 


servations, in fact, indicate that this relativistic expansion occurs along a collimated jet, rather than isotropically 
( [Gehrels et al.l 2009 ). 

If the photons escaped when the high optical depth had decreased to below unity one would have expected 
both a low radiative efficiency (most energy is bulk kinetic) and a quasi-thermal spectrum. Both issues, however, 
are avoided if the radiation arises in shocks, which are expected to occur in the outflow, e.g. |Meszaros| ( j2006| . 
Shocks achieve the reconversion of kinetic energy into random particle energy, and Fermi acceleration in shocks 
leads to relativistic power law particle distributions which in the magnetic fields amplified in the shocks lead to 
synchrotron and inverse Compton non-thermal spectra, when the shock occur beyond the scattering photosphere 
( [Rees & M eszaros, 1992||Rees~& Meszaros, 1994). In addition, a quasi-thermal contribution, possibly with a 
non-thermal tail, can be expected from the photosphere, which however to be important requires internal dissipa¬ 
tion. One expects internal shocks in the outflow when shells ejected with some variability timescale t v with higher 


F catch up with slower F shells around a radius r*- 5 cx 2 ct v r 2 ~ 10 i2 ^ 


IK —3^2 


cm. An external (or termination) 


shock is expected when the jet starts to be decelerated by a swept up external matter of density no (either ISM or 
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a progenitor wind), at a radius r es ~ ( 3Eo/4nnom p c 2 ) l ^r~ 2 / 3 ^3x 10 l6 (E 5 i/no) l / 3 F^ 3 ^ 2 cm, see Figj^] The 
internal shocks (which if magnetic fields are tangled could be a magnetic dissipation region) occur closer in, and 
are highly variable and prominent at y-ray energies (Rees & Meszaros} 1994| ). The external shocks occur further 
out, producing, as the jet is increasingly decelerated, an increasingly softer and longer duration afterglow spec¬ 
trum, from X-rays through optical to radio ( Meszaros & Rees[|1997| . A qualitatively similar spectral-temporal 
behavior is expected for both long and short bursts. 


Fig. 2 
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Such systems, consisting of as much as three different emission zones, the external shock furthermore con¬ 
sisting of both a forward shock and a reverse shock propagating onto the ejecta, can give rise to a variety of 
non-thermal spectra. Upscattering by electrons in the various regions of photons arising in the same region or 
inner regions can give rise to multiple radiation components, e.g. Meszaros] ( |2006| ). 


3 Observational progress through Swift and Fermi 


The observations from the Swift and Fermi satellites, e.g. Gehrels & Meszaros (2012); Omod eTet al.| < |2Q 13 1, 
have led to a number of interesting results. Swift was able for the first time to follow GRB afterglow light curves 
in X-ray and UV/O starting less than a minute after the y-ray trigger. This led, for instance, to the discovery 
in many afterglows of sudden drop of the X-ray light curves at the end of the prompt y-ray emission, followed 
by a shallow decay lasting as much as 500 s in some cases, after which the canonical power law fireball decay 
took over (Fig. [3] left). This behavior appears both in long and short GRBs, as expected from afterglows. The 
sudden drop may be due to the high latitude emission of the jet, while the shallow decay may be due to continued 
outflow, with occasional X-ray flares (as seen in the figure) due to late internal shocks, e.g. jZhang et al. ( 2006 ). 

Swift also made possible an explosive growth in the number of ground-based afterglow localizations with 
redshift determinations (see Fig. [3] right), a smaller number of LGRBs having already been detected earlier by 
Beppo-SAX and HETE-2. Swift, in addition, localized for the first time the host galaxies of short GRBs, which 
are detectable at redshifts typically lower then LGRBs, and confirmed that the hosts of SGRBs include many 


old population (elliptical) galaxies, compatible with an origin in compact double degenerate binaries, e.g. Nakar 
( |2007[ ); |Gehrels et ak| ( [2009| . A significant fraction of short bursts, after their (previously known) initial short hard 
gamma-ray emission, also showed a previously unseen longer (^ 100 s) softer tail in their light curve. Another 
surprise was the first detection of a shock breakout in a supernova accompanied by a low luminosity soft GRB or 
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Fig. 3 Left: schematic X-ray afterglow from Swift-XRT observations (Zhang et all 2006 ). Right redshift distribution of long and 
short GRB (after E. Berger, 2014). 


X-ray flash, e.g. Campana et al.| ( |2006] ); Waxman et al. (2007]). The range of redshifts detected for LGRBs was 
dramatically extended, the highest spectroscopic redshift so far being z = 8.2, with a record photometric redshift 
for GRB090429B at z — 9.4, e.g. [Meszaros & Gehrels| ( |201^| ; |Salvaterra| ( |2015| , see Fig. [3] right. 

Fermi, with its GeV range LAT and MeV range GBM instruments, showed that many, if not most, burst 
detected with the LAT show that the GeV radiation starts to arrive seconds after the MeV trigger, an example 
being shown in Fig. [3] On the other hand, the GeV emission sometimes lasts hundreds of seconds longer than the 
MeV, suggesting association with the afterglow. Various purely leptonic astrophysical mechanisms may give rise 
to this behavior, e.g., |Kumar & Z hang (2014). In addition, this behavior can be used also to provide interesting 
constraints on quantum gravity theories, e.g. Abdo & the Fermi Collaboration (2009b). 



Fig. 4 Light curves of 
GRB080916C show¬ 
ing the GBM (top 
two curves) and 
the LAT (bottom 
three curves) energy 
ran ges |Abdo & the 
IFermi Collaboration 

(2009a). 


Another interesting although not unexpected feature is the presence, in many LAT bursts, of a second spectral 
component at high (GeV) energies, whose slope is harder than the MeV classical broken power law (Band) spec¬ 
tral shape. More interesting was the fact that in some LAT bursts a second separate component is not statistically 
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significant, even though after a few seconds the high energy slope of the Band function (which initially is soft 
and does not extend to GeV) hardens and starts to extend into the GeV range. However, in a larger number of 
LAT bursts a second component with a harder slope does appear, with some delay relative to the softer MeV 
Band function, see Fig. 5. Some theoretical interpretations of these features are discussed in the next section. 


Fig. 5 Spectra of 
GRB090926A from 
Fermi at four different 
time intervals, a= [0.0- 
3.3s], b= [3.3-9.7s], c= 
[9.7-10.5s], d= [10.5- 
21.6s] |Ackermann & 
Ithe Fermi collaboration 

(20 TT ). 
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More recently, a new class of ultra-long GRBs is being recognized. Some of them appear not be connected to 
“usual” GRBs but appear rather to be tidal disruption events (TDEs) of a normal star by a massive black hole at 
the center of a galaxy. The most prominent example of such a TDE detected in y-rays by Swift is the long-lasting 
source Sw J1644+57 (Burrows et al., 201 1) , where the tidal disruption seems to produce a relativistic jet leading 
to y-ray, X-ray, optical and radio emission somewhat similar to a scaled-up and stretched GRB. Other TDEs 

( |2015) , leading to an increased interest in TDE model 


have been detected in X-rays or optical, e.g. 


calculations, e.g. Liu et al. 


|et al.| ( |2015| ); |Mimica et al. 


20121; [Stone et a] 


12015);|Piran et al. 


Komossa 

Mpr 


(2015a 


b). 


Coughlin & Begelman| i |2014| ; |Kelley et al.| i |2014[ >; |Liu 


However some ultra-long GRBs (loosely defined as burst-like objects lasting upwards of 1000 s in gamma- 
rays) may be due to a longer version of stellar events similar to classical GRBs, e.g. Gruber et al. (2011); Stratta 


et al.| ( 20l3] ); Levan et ak| ( |2014) . One way in which such ultra-long gamma-ray bursts (ULGRBs) could arise 


is if the progenitor is larger than the usual Wolf-Rayet progenitor (WRs lost through winds their outer envelope 
and have radii ~ 10 10 — 10 11 cm so jet breakout plus accretion times accord with the classical LGRB durations, 
e.g. |Woosley & MacLadyen (1999)). Instead, ULGRBs might arise in blue supergiants, which have an extended 


(2013); Zhang et al. 

(2014); Evans et al. (2014); Boer et al. (2015); Suwa & Ioka (2011 ); f Nagakura et al. (2012); 

Yoon et al. ( 

2015 

>; 

Gao & Meszaros (2015). Early jet breakout calculations from blue supergiant progenitors 

were discussed in 
were estimated in 

M 

[eszaros & Rees (2001), while ultra-long duration X-ray and O/IR signatures of Pop. Ill stars 

, e.g. Komissarov & Barkov (2010); Meszaros & Rees (2010|); Toma et al. (2011), and the 


possibility of testing for an extended envelope or a Pop III progenitor through a high energy neutrino precursor 
was considered in Razzaque et al. (2003); Schneide r et al.| ( [2002] ); |Gao et al.| ( |201l] ). 


4 Some recent theoretical developments 

The afterglow model based on an external forward and reverse shock is, on the whole, able to explain the basic 
features of the emission in the latter phases. The internal shock model, on the other hand, which is widely used 
by observers as the work-horse tool for interpreting the data, in its simple form has two problems: one is a low 
radiative efficiency, due to the low mechanical shock efficiency, and the other its asymptotic theoretical low 
energy spectrum, which for some bursts is softer than observed. Two main avenues have been explored in the 
past decade to resolve this issue. 






















































































































6 


Meszaros et al. 


Old paradigm: internal + external shock 




< 2005 


Fig. 6 One version of the evolving 
paradigm for GRB emission, where the 
prompt MeV emission is attributed to 
a dissipative photosphere where non- 
thermal processes modify the thermal 
spectrum. 


>2005 


One way to address this is to consider, instead of the usual non-dissipative (and hence inefficient) photosphere 
plus internal shock, a new paradigm consisting of a photosphere which has a hard low energy slope, with nonther- 
mal high energy contributions from Comptonization and/or additional internal shocks ( [Meszaros & Rees 2000). 


The photosphere can be radiatively efficient if it has internal dissipation, e.g. from sub-photospheric shocks (Rees 


& Meszaros |2005( ), and there would still be also an external shock as well, see Fig. [4] Numerical simulations 
(Pe’er et al., 2006| ) including sub-photospheric shocks followed by synchrotron and inverse Compton scattering, 
i.e. a leptonic mechanism, confirm both a high efficiency and reasonable low and high energy slopes as well as 
break energies. 

Another source of sub-photospheric dissipation is inevitable, of the outflow baryonic load contains both pro¬ 
tons and neutrons, which should be the case. This is because the pn fluid, while initially expanding together well 
coupled by elastic collisions, eventually decouple, the neutrons lagging behind, which leads to inelastic collisions 
resulting in pions, whose decay produces both neutrinos and relativistic positrons (Derishev et al., 1999; Bahcall 


|& Meszaros) |2000| ). This also results in efficient dissipation, the non-thermal positrons heating the plasma and 
leading to escaping photon spectra which again are reasonable (Beloboro dov||2010| ), see Fig. 7. 



Fig. 7 Dissipative photosphere 
heated by np collisions leading to 
a non-thermal “Band” spectrum ( [Be-| 
|loborodov] 12010) . The high energy 
slope extends into the GeV band thanks 
to energy injection from the non- 
thermal decay positrons. In this case 
it reproduces a single Band function, 
but adding magnetic fields can add a 
second high energy component with a 
harder slope ( |Vurm et aL]|2013) . 
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A more robust way to obtain a second high energy component in such dissipative photospheric models is 
through the up-Comptonization of the inner photospheric spectrum by the external forward and reverse shock 
(|Veres & Meszaros||2012), see Fig. [8] 

GRB 09051OA 




Fig. 8 Left: schematic dissipative photosphere and upscattering by external shock. Right: baryonically dominated dissipative photo¬ 
sphere plus external shock upscattering model leading to a two-component Band plus GeV spectrum ( |Veres & Meszaros||2012) . 


The other main way to address the efficiency and spectral issues of internal shocks is by including the accel¬ 
eration of hadrons in the same internal shocks, which is ignored in the classical internal shock picture, but which 
should be essentially unavoidable, if the outflow is baryonic. Hadron acceleration in internals shocks was first 
considered by Waxman & Bahcall ( |1997| ) to predict TeV neutrinos from GRBs, and the associated high energy 
y-rays were considered by various authors, e.g. Bottcher & Dermer (1998); Asano et al. ( 2009b). The high energy 
components are produced by the synchrotron and inverse Compton of the cascade secondaries, see Fig. [9] left. 
The delay associated with the acceleration and cascade development can also possibly explain the delay of the 
GeV emission relative to the MeV radiation ( [Asano et aL||2009b]|Razzaque||2010tpVsano & Meszaros|[2013| ). 


Simple hadronic internal shock cascade model 




Fig. 9 Left: schematic spectrum internal shock with electron plus hadron acceleration and secondary cascades. Right: numerical 
calculation of hadronic internal shock leptonic plus cascade secondary spectrum (Murase et al. 2072) 
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The hadronic secondaries can reproduce both the high energy second component, and also an additional low 
energy component which appears associated with the high energy component. This is illustrated as the leftmost 
component in Fig. [9]left, a hint of which is seen in the low energy portion of the data of Fig. 5, and more clearly 
in GRB090510, a calculation of which is given in Asano et al.] f2009a] ). This component could also give rise to the 
bright early optical flash seen in some bursts. An interesting possibility is that re-acceleration of the secondaries 
via Fermi 2nd order mechanism in the turbulence behind the internal shocks can lead to a self-consistent Band 
spectrum at MeV energies, see Fig. [9] (right), which thus may not need to be due to the photosphere. The radiative 
efficiency in this case is also high, as well as the spectrum being reasonable, without an unduly high proton 
luminosity. 


5 Multi-messenger approaches 


The multi-messenger signatures of GRB are an important direction being actively pursued through multiple 
channels, both photonic and non-photonic. The latter, e.g. in the form of neutrinos, cosmic rays and gravitational 
waves, are being pursued by multiple observatories, such as IceCube, Auger, Telescope Array, LIGO and VIRGO, 
in addition to VHE y-ray facilities such as HAWC, VERITAS, HESS and others, complementing Fermi, Swift 
and other existing instruments. An ambitious effort to coordinate the sub-threshold signals from various such 
facilities is the AMON network (http://amon.gravity.psu.edu, see [Smith et al. ( 2013| ). An interesting question is, 
if GRBs accelerate hadrons, can they explain at least partly the Auger CR and IceCube neutrino observations, or 
do they satisfy the constraints imposed by these observations? And will short GRBs, if they are indeed merging 
neutrons star binaries, be detected in the near future by gravitational wave observatories? 

With IceCube and Auger in full operation, the first question is being already addressed. Auger has measured 
the ultra-high energy cosmic ray (UHECR) spectrum to great precision (The Pierre Au ger Collaboratio n et al. 


2013| ), but since CRs are deflected by intergalactic and galactic magnetic fields, one can only test whether GRBs 


can produce the appropriate flux within the GZK radius, as suggeste d by|Waxman (1995]). Since the GRB UHE- 
CRs would produce via py TeV-PeV neutrinos (Waxman & Bahcall, 1997), these should not exceed the IceCube 
limits. The first limits were set by the IceCube initial 40+59 string array ( Abbasi et al.||2012~] ), which compared 
neutrino observations at the locations and time windows corresponding to 300 GRBs detected electromagneti- 
cally by satellites, 183 GRBs with the 50 string array and 177 GRBs with the 40 string array. The comparison 
used the simplified internal shock photo-hadronic model of [Waxman & Bahcall] ( [1997] ) to predict an expected 
neutrino flux, resulting in an over-prediction of this model compared to the data by a factor ~ 5. A more generic 
model independent test, where the detected Auger cosmic rays are attributed to the GRBs was also negative. 
However, more accurate descriptions of the photo-hadronic interactions as well as inclusion of astrophysical 
model uncertainties ( |He et al.| |2012| |Li| 2012| |Hummer et al.| 2012| ) showed that discrepancies between this 
model and IceCube might only be expected after several years of observations with the full array. These tests 
also used stationary radius internal shocks, allowing at most for range of radii (through t v ). 

However, the internal shocks, after forming, continue to move out, and doing so expand, so the time depen¬ 
dence of the varying comoving quantities must be taken into account. Such a time-dependent calculation was 
performed by [Asano & Meszaros] ( |2014| ), who took into account also the observed range of t v as well as the 
shock expansion and the GRB luminosity function ( Wanderman & Piran, |2010| ) to obtain the expected flux of 


y-rays, neutrinos and neutrons escaping from the GRBs - see Fig. 10 (right). 


These sources were used then to calculate the expected diffuse fluxes, using the luminosity function of Fig. 
10] (left). These predicted diffuse fluxes are compared in Fig. [Tl] to the IceCube 40+59 limits ( Ab basi et al.| 


2012) and the diffuse Auger UHECR flux ( [The Pierre Auger Collaboration et al. 2013| . Two prescriptions were 


used for the UHECR escape, one (Fig. |TT] left) assuming that they escape as neutrons due to py interactions and 
then reconvert to protons once outside the source; the other extreme being simple free escape as protons (Fig. 
[IT] right). It can be seen that such GRB time-dependent internal shocks can explain at least part of the highest 
energy UHECR in the GZK range, without violating the 40+50 string IceCube neutrino limits. 
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Asano & Meszaros,ApJ 785:54 (2014) 


ef(e) [erg] 5t=0.1 s 



Wanderman & 


Fig. 10 Left: observed distribution of variability times t v (left top) and GRB luminosity function (left, bottom, from 
Piran (2010). Right: Time-dependent hadronic internal shock Monte Carlo cascade spectral calculation for a distribution of variability 
times and observed luminosity function, showing photon, neutron and neutrino spectrum ( [Asano & Meszaros||2014) 



Fig. 11 Left: the CR (black) and neutrino (red) diffuse intensities for the neutron conversion model, compared to the Icecube 40+59 
(Abbasi et al. 2012) and the Auger data (The Pierre Auger Collaboration et al. 2013). Gray shaded area for CRs indicates the 
uncertainty in the local GRB rate. The thin dash-dotted line is the CR spectrum without the effects of photomeson production and 
Bethe-Heitler pair production. The cosmogenic neutrino spectrum produced via the GZK process is also shown as the thin red line. 
Right: same as left, but for the sudden release model. From |Asano & Meszaros| ( |20T4| 


As mentioned earlier in Q one of the new approaches to describing the prompt GRB y-ray emission ascribes 
it to a dissipative photosphere, with a diminished role for the internal shocks (e.g. if the outflow is magnetically 
dominated, internal shocks might be prevented). Thus, it is necessary to consider also the possible hadron accel- 


eration and neutrino production in the GRB photospheres (Murase 2008 

Wang & Dai 2009; Gao et al. 201^ 

Asano & Meszaros, 2013; Zhang & Kumar, 2013; Murase et al. 2013b;; 

Bartos et al. 

2013). In addition, modi- 

fled internal shocks might also include magnetic dissipation, e.g. the ICMART model ( 

Zhang & Yan 2011), and 


both these and photospheric models predict qualitatively different neutrino fluxes than classical internal shocks. 
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A recent IceCube analysis, using 4 years of accumulated data including now the full array, tests such models 
( [IceCube Collaboration et akj |2014| ), concluding that less than 1% of the electromagnetically detected GRBs 
(modeled as either classical internal shocks, ICMART shocks or photospheres) can be contributing to the diffuse 
IceCube neutrino flux. 

An exciting development is the discovery by IceCube of a diffuse neutrino flux which appears to be of astro- 
physical origin (Aartse n et al. 2013[ IceCube Colla boration 2013|), with limited (or no) contribution from our 


own galaxy ( |Fox et ~ 2013 11 Ahlers & Murase 2014[ Anchordoqui et al.||2014a||Ahlers et al.||20T5] ). This flux 
has a spectrum extending from about a few PeV down to ^ 10 TeV, with a slope which (if unbroken) appears to 
be ~ —2.3 to —2.5. As mentioned above, less than 1% of these can be coming from “classical”, i.e. normal high 
luminosity GRBs detected by Swift or Fermi ( [IceCube Collaboration et al.[|2014| ). In fact, looking at Fig. |TT] it 
is clear that such GRBs, while not violating the IceCube observational limits, in fact should not be producing 
anything near the observed neutrino flux in the TeV-PeV range. The question then is, what are the sources of 
these PeV neutrinos? If they are connected to GRBs at all, it might perhaps be to “hidden” GRBs, which are 
electromagnetically weak and do not trigger Swift or Fermi ([Murase & Ioka||2013 Liu & Wang||2013| ); these 
could be either so-called choked-jet GRBs (Mesz aros & Waxman 2001), whose jets never emerged from the 
progenitor star, or they might be (electromagnetically) low-luminosity GRBS, or LLGRBs, e.g. (Mu rase et al.| 
2006; |Gupta & Zhang||2007[|Murase et~aL[|2008 ). 

Other, non-GRB related possible sources of the diffuse TeV-PeV neutrinos could be galaxy cluster accretion 
shocks and/or starburst galaxies, e.g. |Loeb & Waxman| ( |2006| ); [Thompso n et al.| ( [2006| ); |Murase et al.| ( |2013a| ); 
He et al.| ( |2013| ); [Anchordoqui et al. ( 2014b ); Tamborra et al.|(2014 ); Chang & Wang ( 2014|); Chang et al. ([2014); 


or they could be due to AGN jets, e.g. Becker Tjus et al. "([2014 ); Murase et a l 
they could be due to galaxy-galaxy collisions ( [Kashiyama & Meszaros 2014 ) 


( 2014| ); KrauB et al. ( 2015| ); or 



Fig. 12 Diffuse neutrino flux per flavor (solid 
black) and y-rays (dash-dot) from HNe and SNe, 
for a diffusion coefficient /)«£/ , in both the 
host galaxy and cluster. The proton spectral in¬ 
dex is r = 2. The black line with white cir¬ 
cles denotes the measured atmospheric neutrino 
flux. The starburst galaxy contributions are in 
red lines, while the normal star-forming galaxy 
contributions are represented by blue lines. The 
contribution from HNe are in solid lines colored 
while those from the SNe are dashed. The solid 
green line denotes the total cluster contribution 
(i.e. HNe and SNe from both types of galaxies). 
Green data points correspond to the Fermi mea¬ 
surements of the extragalactic diffuse y-ray back¬ 
ground ( [Ackermann et al.| |2015| . Black points 
correspond to the IceCube measurements of as- 
trophysical neutrinos (Aartsen et al. 2015), note 
that two of the low energy data points are within 
the gray lines of the error bars of the atmospheric 
flux. From |Senno et al.| ( [2015) 


On the other hand, a possible UHECR and TeV-PeV neutrino source closely related to GRBs are hypernovae, 
or jet-driven supernovae ( Wang et al.||2007} Budnik et al.| |2008 ), and these are in fact largely concentrated in 
star-forming or star-bursting galaxies. While classical long GRBs are often associated with hypernovae, there 
are many more (90-95%) of hypernovae which do not show signs of a GRB, i.e. non-GRB hypernovae are 
much numerous than GRBs. Hypernovae have also been investigated as possible sources of the IceCube diffuse 
astrophysical TeV-PeV neutrino flux (He et al.| |2013| |Murase et al.| |2013a[ |Liu et al] |2014| ). The challenge, 
with these various sources, is to fit the IceCube diffuse spectral flux, arising from pp or p y collisions leading to 
7T + , while at the same time not violating the observed Fermi isotropic diffuse y-ray flux (Abdo & Fermi LAT 
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Collaboration, 2010; Ackermann et al.||2015] ). These constraints can be reasonably met when taking into account 
both the hypernova and normal supernova populations in the observed relative ratios, as seen in Fig. 12, ( Senno 
et al. 2015[ ). The protons accelerated in the hypernovae (HNe) and supemovae (SNe) propagate in the host galaxy 
and in the host cluster with a diffusion coefficient D(e) oc e 1 / 3 or oc r 1 / 2 , where £ is proton energy, and make 
neutrinos and gamma-rays mainly by pp collisions, with negligible py contributions. As can be seen in Fig. 12, 
the IceCube flux may be satisfied (except for one point a few sigma away) if the entire Fermi flux is explained this 
way. The problem is that AGNs already appears to be responsible for a large fraction of the diffuse Fermi y-ray 
flux (Abazajian et al.| |2011[ ), hence not all the IceCube flux may be ascribed to HNe and SNe. Future work needs 


to clarify the roles of the various possible sources. The new all-sky y-ray survey instrument HAWC (Westerhoff 


2012 


2014) will undoubtedly play a large role in the investigation of GRBs (Abeysekara & the HAWC collaboration 


Taboada & Gilmore} |2014[ Abeysek ara et al.| |2015| ) and other extragalactic as well as galactic sources of 


cosmic rays and neutrinos. 
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